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Abstract 

We observe that probing certain classical field theories by external sources uncovers the un- 
O ! derlying renormalization group structure, including the phenomenon of dimensional trans- 

mutation, at purely-classical level. We perform this study on an example of X(f)^ theory and 
unravel asymptotic freedom and triviality for negative and positives signs of A respectively. 
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H I We derive exact classical (3 function equation. Solving this equation we find that an iso- 

lated source has an infinite energy and therefore cannot exist as an asymptotic state. On 
the other hand a dipole, built out of two opposite charges, has finite positive energy. At 
large separation the interaction potential between these two charges grows indefinitely as a 
distance in power one third. 



I. INTRODUCTION 

The discovery of asymptotic freedom in QCD [l| opened a new era in particle physics. 
Besides its direct relevance for understanding the nature of strong interactions, it introduced 
a concept of dimensional transmutation or equivalently a dynamical generation of scale 
in a seemingly scale-free theory. Both phenomena, asymptotic freedom and dimensional 
transmutation, are usually perceived as intrinsically quantum phenomena, as they are both 
deeply rooted in the renormalization group properties of the quantum theory. The question 
we would like to address in this paper is, how much of these phenomena is captured by 
classical physics? This is a fully legitimate question, since usually quantum effects have 
classical precursors, which sometimes appear in the form of uncontrollable growth of the 
classical fields. 

In order to illustrate our ideas, we will consider a simple example, namely, scalar theory 
with negative A0^, which is known to be renormalizable and has a negative /3 function 
|2|, but has an unbounded from below potential. However, our classical renormalization 
group treatment delivers a natural prescription which allows to self-consistently work with 
this theory in the presence of external sources, and isolate our findings from the issue of 
potential instability in a pure A0^ theory. 

In the presence of the external sources, the requirement of independence of physical 
observable from the regulator scale of sources implies the running of the effective coupling. 
This running is the main reason behind the whole renormalization group structure and 
the subsequent classical dimensional transmutation. This scale dependence is the key to 
why the renormalization group results can be safely disentangled from the instability issues. 
Simply speaking, because of the emerging scale-dependence we always perform calculations 
on time-scales shorter than would-be instability time in a sourceless theory. 

Surprisingly, by probing X(j)^ theory by the large external source, we uncover the whole 
built-in RG structure already at the classical level, with fuUy-fiedged counterparts of asymp- 
totic freedom as well as dimensional transmutation phenomena, in which an analog of QCD- 
scale appears as a result of classical RG invariance. We derive the exact classical (3 function 
equation from which we extract non perturbative information about the infrared region. 
Solving this equation in strong coupling regime we find that the energy of the isolated ex- 
ternal charge is infinite and positive and hence it cannot exist as a free asymptotic state. 



Moreover, considering a dipole, build out of two charges with opposite signs, we find that 
its energy is positive and finite. When charges in this dipole are separated their interaction 
potential grows indefinitely as distance in power one third, thus confining the charges. These 
findings indicate that there may exist a classical counterpart of confinement. 

II. CLASSICAL SOLUTION 

After the discovery of asymptotic freedom in the non-Abelian Yang Mills theories, the 
physics underlying the negative sign of the /3 function was understood as an ant i- screening 
effect due to self-interactions of the gauge fields. From the classical point of view we can try 
to understand this anti-screening considering how self-interaction modifies the field created 
by an external point-like source Q at large distances. This modification compared to the case 
of free fields can be used to define the effective charge Qeffij-) at distance r or, equivalently, 
the running coupling constant a{r). This can be done without invoking quantum theory 
and the result will only depend on the particular classical features of self-interactions. For 
the large external charge one can expect that the classical contribution to anti-screening will 
dominate over the one due to the vacuum polarization effects. 

In this section we will solve perturbatively the classical equations of motion for a given 
external charge and show how the anti-screening effect (the growth of Qeff{i^) with r) can 
naturally be achieved. 

A. Anti-screening 

Let us consider Ao0*^ theory with negative Aq 

S = J (^d^(pd^<P - ^Ao0' + AtiqA d'x, (1) 

where the signature is taken to be +, , and Q is the external charge. In the case of 

a point-like charge the field equation for the static spherically symmetric field (r) reduces 
to 

4#^'#V^o0^ = -4^^^W- (2) 

r"^ ar \ dr J 

If one neglects the nonlinear term in this equation then its solution is 

00 = ^. (3) 

r 



Notice that for Q ^ 1 the amphtude of the classical scalar field on scales r is much larger 
than the typical amplitude of the quantum fluctuations, which is of order 1/r. If the coupling 
constant Aq is small enough, i.e. XqQ^ <^ 1, then the corrections to solution ([2]) due to the 
self-interaction Ao0^ can be treated perturbatively. The leading order correction to the 
solution 00 can be obtained by solving equation 

^^(y-^)=-4.QH-)^^.%, (4) 

r^ ar \ dr J v^ 

where we have substituted 0o in the nonlinear term. The last term in this equation can be 
treated as the contribution to the effective charge induced by the nonlinear self-interaction. 
As it was noticed above, for Q ^ 1, the vacuum polarization contribution to the induced 
charge is much smaller than the classical contribution and therefore can be neglected. As 
one can easily see from (jl]) the induced charge for positive Aq has a sign which is opposite 
to the sign of the source and the nonlinear interaction leads to screening. For negative Aq 
the charges have the same sign and we have an anti-screening effect similar to the one of 
the non-Abelian gauge theories. Since we are mainly interested in asymptotic freedom and 
confinement we will consider only the case of negative Aq. 

B. The perturbative expansion 

Let us look for the solution to equation ([2]) in the following form 

Integrating equation ([2]) and substituting this ansatz we can rewrite the equation for the 
scalar field as 

/ (r) = 1 + aor ^ (^/^ ^^-" j -, ' N (a,) , (6) 

where 

ao = -AoQ' > 0, (7) 

is the effective coupling constant and we have introduced the ultraviolet cutoff scale tq to 
regularize the integral, which otherwise would diverge. The function A^ (oq) ? which depends 
only on oq, is fixed by the normalization condition: / (ro) = 1. It is clear that in the limit 
ao — )■ it must vanish and therefore in the absence of self-interaction, the solution ([5]) 



with / (r) = 1 exactly satisfies equation ([2]). As it follows from ([5]), the function / defines 
the anti-screened effective charge Qeff (r) = Qf (r) or, equivalently, the running coupling 
constant 

"e// (r) = aof (r) (8) 

Assuming that ao ^ 1 we can solve the integral equation ([7]) by iterations in powers of oq. 
With this purpose it is convenient to rewrite it as 

f{x) = l + aoe" [ f W) dx" e-""' dx' - N (ao) , (9) 

where we have introduced x = In {r/ro) instead of r. Substituting / (x) = 1 into the right 
hand side of equation ([9]) and taking into account that A^ (ao) = ao + O (oq) "^^ fi'^d 

f{x) = l + aQX + [al) . (10) 

Next we take this solution, substitute it again in ^ and take A^ (ao) = ao + Sag + O (oq) . 
Keeping only the terms up to second order in a^ leads to 



Qx^ + 3x^+0 [al) . (11) 



/ (x) = 1 + a^x + Oq ( 2^ 

This procedure can be iterated giving us at each step the next order term in a^. The result 
up to order a^ is 



/ (x) =1 + aox + ali-x^ + d,x\+ali -x^ + 12x^ + 24x J 
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+ al I —x^ + —x^ + ^^x^ + 285x 



+ al ( —x^ + 93x^ + — ^x^ + 2142x2 + 4284x 



+ al i—x^ + ^^x^ + ^^x^ + 10 521x3 + 37 989x2 ^ 75 g7g^^ _^ ^ / 7n _ 
" V 16 40 4 y V u/ 

(12) 
The function A^ (ao) to the same order in perturbations should be taken to be 

N (ao) = tto + 3^0 + 24^0 + 285a^ + 4284a^ + 75 978a[; + O [al) . (13) 

The effective running coupling as a function of distance can be written then as perturbative 

series in powers of ao 

00 
«e// (x) = aeff (r) = ao/^ (r) = ^ a^"^^^„ (x) , (14) 

n=0 



where 

go (x) = 1, gi (x) = 2x, g2 (x) = Ax^ + 6x, g^ {x) = 8x^ + 30x^ + 48a;, 

c/4 (x) = IQx^ + 104x^ + 342x2 + 570x, 

^5 (x) = 32x^ + 308x^ + 1572x^ + 4998x2 + 8568x, 

4176 
ge (x) = 64x^ + ^— x^ + 5880x^ + 27 612x^ + 86 832x2 + 151 956x, (15) 



etc. The calculation of gn (x) is straightforward and we did it until gio (x) . However, to 
simplify the formulae we present here the result only up to gg (x) . Note that the running 
coupling constant depends on r only logarithmically with the coefficients gn power series 
of X = In {r /tq) with the highest power n. Moreover, the series flT^ can be rearranged 
and partially resummed. In particular, collecting together leading powers of logarithms, 
next-to-leading and next-to-next leading powers we get 



a. 



eff (x) = ao [1 + 2x + 4^2 + 8x^ + 16x^ + 325^ + O (x^)] 



+ «o 
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6x + 30x2 ^ ^Q4~3 ^ 308^4 ^ HZ^^s ^ HZ^^e ^ q (^~7^ 
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48x + 342x^ + 15725'' + 5880x^ + 



97248.. 



-X 



+ 59 248x'' + O (x^) 



+ 0K..) 



(16) 



where x = a^x. In the second and third brackets we have also included higher order terms 
compared to (TT^ . The series in the bracket can be resummed. In particular, it is obvious 
that 



1 + 2x + 4x2 + 8x^ + 16x^ + 32x'' + ... 



1 



l-2x' 



Much less obvious are the following results 
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(18) 
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48x + 342x^ + 1572x^ + 5880F + 



^.3 



97248 



^5 



x° + 59 248F... 



9 (In (1 - 2x)y - 9 In (1 - 2x) + 30x 



(19) 



(l-2x)' 

which the reader can verify just expanding the appropriate expressions in powers of x. 
One may wonder how did we manage to resum these last two series? The answer to this 



question is in the next section where we uncover the renormahzation group structure of our 
entirely classical theory and derive the /3 function which generates the resummation of the 
perturbative expansion to the appropriate powers of cto- 

III. RENORMALIZATION GROUP AND ASYMPTOTIC FREEDOM 



In the Wilsonian approach [3j the renormahzation group sets how the couplings of the 
quantum theory should change under re-scalings of the UV cutoff . The equations governing 
this dependence are known as the renormahzation group equations. This general notion of 
renormahzation group can be extended to the classical field theory with external point-like 
sources in the following sense. Let us introduce an UV cutoff tq setting the way we smear the 
source. The classical field created by such source will generically depends on the regulator 
To and the self-coupling A of the theory. One can ask under which circumstances we can 
require that the classical theory must be invariant under re-scaling of tq. This is possible 
only if the corresponding classical theory incorporates the renormahzation group structure. 
In this case the dependence of the coupling on the smearing cutoff A(ro) also captures the 
screening and ant i- screening effects. Moreover, using the effective running coupling we can 
associate with an external source, a physical length scale Re by the standard procedure of 
dimensional transmutation. 

The reason why the classical solution captures the renormahzation group structure is 
easy to understand. Any regularization scheme in quantum field theory give rise to log- 
arithmic contributions which even in a scale invariant theory lead to anomalous scaling. 
These logarithmic contributions are of the type log (p/A) with A the UV cutoff. Since the 
divergent contribution log A is absorbed by renormahzation, we are free to choose the scale 
p at which the logarithmic contribution to the self-energy vanishes. As a consequence the 
scaling of A should be accompanied by finite renormalizations (the RG transformations) of 
the coupling constants. In the classical theory under consideration we have found the same 
type of logarithmic contributions to the field created by the external source. In this case the 
role of cutoff A is played by the smearing scale tq. One can renormalize the classical theory 
by subtracting the log Tq contributions as it is done in quantum field theory. However, if we 
want physics to be independent on the method of removing this infinity, we need to change 
the couplings, exactly as it is done in quantum field theory. Both renormahzation group 



structures, the classical and the quantum ones, are structurally identical for Xcp^ theory 
because both have the same type of parent logarithmic contributions. 

A. Perturbative /? function 

Once we have obtained the classical expression (fT^ for the effective coupling constant 
Oeff (r) we can check whether taking the bare coupling constant oq as a function of r^ we 
can make a^jf (r) independent of tq. As we have said, this is possible only in the theories 
with associated renormalization group structure, which in turn imposes rather severe con- 
ditions on the functions Qn (x) in the perturbative expansion flT^ . Let us first derive these 
conditions, which do not depend on the origin (classical or quantum) of the renormalization 
group, and then verify whether they are satisfied by the functions in flT5|) . 

On general grounds the expansion of the dimensionless running coupling constant Og// (r) 
in powers of ao = a (^o) , normalized at r = ro, can be written as 

ae// (r) = a (ro) + a^ (ro) g, ( -] + ... = V «"+' (ro) gn ( -] , (20) 

where we use the spatial scale r instead of the usually used energy scale k ~ 1/r. It is clear 
that go (r/ro) = 1 and since a{r) = a (ro) at r = ro, we have 

9n{l) = 0, (21) 

for n > 1. Invar iance under changes of the ultraviolet regulator ro, implies 

which in turn imposes severe restrictions on (?„ (r/ro) . Taking the derivative and rearranging 
the terms in (!22|) leads to 

oo 



E 9k+i (^) " (^o) 
«' (^o) -.^ , (23) 



da{ro) 2/ N k=o 



k=0 

where a; = In (r/ro) and prime denotes the derivative with respect to x. The ratio of sums in 
the right hand side of ( 123|) should not depend on x because the left hand side of this equation 



is X- independent. Therefore setting x = (which corresponds to r = ro) and taking into 
account f l2T|) we find that 

oo 

-^^ = E aUi (0) «' (ro) , (24) 

Eik + l) gk (x) a'' (ro) fc=o 

fc=0 

from where it follows that the functions Qk {x) should satisfy the following recursion relations: 

dgn+i (x) 



J2ik + l)9'n+i-k{0)9k{x). (25) 



dx 

k=0 

Note that only if these conditions are satisfied then there exists a function a (ro) for which 
the sum in the right hand side of (120|) does not depend on tq. Nicely enough the unambiguous 
solution of these recursion relations with "initial conditions" ( 12T|) is given by 

n 

9nix) = J2ckx'', (26) 

fc=l 

where c^ are completely determined by the numerical values of g[ (0) , (^2 (0) , •• which in 
principle can be arbitrary. For instance, the coefficient in front of the leading logarithm 
x^ = In" (r/ro) in gn is equal to c„ = {g[ (0))" . 

At this point it is quite rewarding to confirm that the set of classical functions ( ITSl) in fact 
satisfies the recursion relations ( l25l) . This can be done by direct calculation to any order in 
perturbation theory (we did it up to gio). Thus, taking oq in (fT4|) to be the function of ro 
we uncover the renormalization group structure of the classical A0^ theory. We would like 
to stress that in distinction from the quantum field theory, where the renormalization group 
is normally checked by direct calculations only to the leading logarithms (and postulated 
otherwise), we verified it also for all subleading logarithms. 

To take the advantage of renormalization group for partial resummation of the perturba- 
tive expansion (fT^ we note that from (l23l) and fl2^ it follows 

da (ro) 



dlnro 



a' 



(ro)E^U(0)«'=(ro). (27) 

fe=0 

The running constant Og// (r) depends on r in the same way that a (ro) depends on rg. 
Hence aejf (x) satisfies the well known Gell-Mann-Low equation [^ 

^^^ = «e// (^) E ^^+1 (0) <f (^) • (28) 

fc=0 



The P function is normally defined as the derivative of Oe// with respect to the logarithm 
of the energy squared. For us it is more convenient to define it as 

13 ^ ^^#^, (29) 

ax 

which (up to factors 4tt due to the choice of charge units) is related to the standard (3st 

function as f3st = — /3/2. According to ( !28l) and ( fTSl) the classical perturbative (3 function is 

equal to 



f3{a) = J2 Pka''^^ = 2a^ + Qa^ + 48a^ + 570 



a^ 



fc=i 



+ 8568a^ + 151956a^ + ..., (30) 

where a = a^fj (x) and /3i = g[ (0) . 

Obviously, we should not expect the numerical coefficients /3j of this classical beta func- 
tion to be identical to the ones derived in the quantum field theory. In the last case /3j 
are determined by the loop contributions and they will depend, beyond two loops, on the 
particular renormalization scheme used to segregate a finite part of the divergent loop inte- 
grals. The classical beta function accounts for the ant i- screening effects due to the classical 
self-interaction. A potential quantum theory check of the numerical coefficients derived 
above will require to work in the presence of large external charge where we have to modify 
the Green functions in order to account for the effect of the external charge. Because for 
Q ^ 1 the quantum fluctuations are subdominant we expect that classical contribution 
dominates. Although the direct check of this expectation is obviously important we will 
not follow that path. Instead we will restrict ourselves to the physical consequences of the 
underlying renormalization group structure of the classical theory. 

B. Partial resummations 

In equation ( 1TB]) we have separately collected the contribution of the leading and sublead- 
ing logarithms to a (x) and presented the result of their resummation. For the subleading 
logarithm the result was derived using Gell-Mann-Low equation. For a ^ 1 we can flrst 
neglect all terms in /3 function besides of the "one loop" contribution. Equation fl28|) then 
reduces to 

^ = 2a^ (x) , (31) 

10 



and its solution, with initial condition a (0) = «o = —^oQ^, is 

/ \ _ Q'o —^qQ /oqn 

" ^"'^ ~ 1 - 2aox ~ 1 + 2AoQ2 In (r/ro) ' ^ ' 

where Aq = A (ro) < 0. It is easy to see that this solution gives us the resummation of the 

leading logarithms in the expansion flT6|) , (fT7|) . 

We can repeat the same analysis keeping in Gell-Mann-Low equation the contribution up 

to two loops, 

^^ = 2a'{x) + Qa\ (33) 

ax 

Integrating this equation with initial condition a (0) = ao? we obtain 

1 3 j^ p + 3a(x) ^ ao \ ^ l-2aox_ ^^^^ 



a{x) \ 1 + Sao a{x)J Oq 

Solving this equation in terms of the perturbative expansion in oq one gets 

Note that the second term agrees with resummation (fTS!) of the next to the leading order 
logarithms. The same is true at three loop order, where the solution to 

da (x) 



dx 
which is 



2a^ (x) + 6a-^ + 48a^ (36) 



a(x) = sf ^ ln(l-2«ox) (37) 



+ 9 (t-^ ) V""^ ^^ ~ ^"0^^ - In (1 - 2aox) + ^aoxj + O (a^) , 

also accounts for the resummation of next-to-next subleading logarithms. In other words 
the solutions to the classical renormalization group equation give us the resummation of 
the perturbative series ( TT^ taking care in every step about next logarithms in gn (x) . It 
is clear that when the running coupling constant becomes of order unity (strong coupling 
regime) all terms in expansion (!37|) are of the same order and the series (j37|) should be 
further resummed. It is not a priori clear whether the singularity in this expansion (Landau 
pole[5J]) will survive after this resummation. We will answer this question in the next section 
using nonperturbative methods. 
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C. Dimensional transmutation and asymptotic freedom 

One important consequence of the renormalization group is dimensional transmutation. 
We can easily understand this phenomenon using the result of the one loop resummation 
of perturbative expansion 

In this expression Aq depends on the regulator tq in such a way that a (r) is ro-independent 
to the corresponding order. Therefore we can define the renormalization group invariant 
physical scale R^ via 

Note that this scale 

i?^ = roe"2M'-o)Q% (40) 

does not depend on the particular value of regulator tq at one loop level. Using this dynam- 
ically generated scale we can rewrite the physical running coupling as 

"«"-^ ''■'«' =2h7(kM- (^') 

The physical meaning of this expression is obvious. Perturbatively the theory can be defined 
in the ultraviolet region corresponding to length scales r ^ R^, where it becomes effectively 
free. Thus, we have found asymptotic freedom in the classical A0^ theory with negative A. 
In the infrared at length scales of order Re the theory becomes strongly coupled and 
non-perturbative. What is the potential meaning of this dynamically generated scale? From 
the point of view of the classical theory the existence of this scale is quite surprising since 
it is independent of the UV regulator. On the top of that Re is a very non-perturbative 
scale. Obviously it is tempting to think of Re as setting the natural confinement scale of 
the theory. A way to check this claim is to derive an exact classical (3 function equation 
and to read off the previous perturbative expansion from the corresponding solution of this 
equation. We address these issues below. 

IV. BEYOND PERTURBATION THEORY AND ASYMPTOTIC BEHAVIOR 

The usual way to address the non-perturbative phenomena within perturbation theory is 
to study the convergence of the perturbative series. In p it was found that the numerical 

12 



coefficients in the perturbative expansion of f3 function asymptotically grow as f3k ~ kipf, 
where /3i is one loop /3 function and k denotes the perturbative order in coupling constant. 
Such behavior sets the limit of perturbation theory and fixes the uncertainty of the computa- 
tions to be of order exp (— l//3ia) . In the theories with asymptotic freedom this uncertainty 
is extremely small in the deep UV region contrary to what happens in the theories with UV 
Landau pole. Normally in quantum field theory is hard to prove this factorial asymptotic 
behavior of the coefficients in /3 function. It can have different origin: either the growth 
of the number of diagrams contributing to a given order in perturbation theory (instanton 
effect) or the contribution of multi-bubble diagrams (renormalons) . 

In the classical approach to the renormalization group, however, there is an opportunity 
to convert the classical equations of motion into exact equation for the /3 function. This 
equation be can used afterwards to check the similarity between classical and quantum renor- 
malization groups. In particular, as we will see, one can use the exact equation to derive the 
asymptotic behavior of the coefficients in the perturbative expansion of f3 function. Inter- 
estingly enough the asymptotic behavior anticipated by the exact classical /3 function agrees 
with the quantum filed theory expectations. In addition this allows us to clarify the origin 
of renormalons as well as the generic form of the non-perturbative uncertainties. Moreover, 
the non-perturbative contributions will be naturally defined in terms of the dynamically 
generated scale Re, as it should be. 

A. Exact classical /3 function 

To derive an exact equation for the classical /3 function we begin with equation for static 
scalar field outside an external source 



Substituting 



5^, (43) 



we can rewrite the equation above as 

/" - /' + aof = 0, (44) 
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where ao = — AqQ^ and prime denotes the derivative with respect to x = In (r/Rc) . Multi- 
plying this equation by ao/ and defining running coupling constant as before 

a (x) = aof (x) , (45) 

we obtain the following second order differential equation for a (x) 

a'2 
a" - -a' + 2a^ = 0. (46) 

2a 

Recalling the definition of /3 function, (] = a', and taking into account that 

a" = 01^4- = ^f^ (47) 

da da 

equation fB6l) reduces to the first order differential equation: 

/5 = 2^2 ^ - -^ - ^ , (48) 



2 \da a ^ 
which determines the exact classical /3 function. 

B. Weak coupling expansion and renormalons 

Let us first use the exact equation for f3 function to reproduce our perturbative results 
above. In order to do that we substitute in (1^81) 



P{a) = J2Pka''+\ (49) 

k=l 

This leads to the following recursion relations for the unknown numerical coefficients (3^ 

I3i = 2, p, = J2[m + M f3k-mf3„, for k>2, (50) 

Using this relations we find 

/3i = 2, /32 = 6, (33 = 48, (3^ = 570, (35 = 8568, /^e = 151956, ... (51) 

in complete agreement with ( !30l) to an arbitrary order in a. Thus we have proven that equa- 
tion (jlSj) yields the exact (3 function which is in complete agreement with the perturbative 
(3 function. 
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Let us now find the asymptotic behavior of the perturbative series. One can easily see that 
for large k the main contribution to the sum in fj50|) comes from the terms with m = k — 1 
and m = 1 and the recursion relation reduces to 

/3fc~(fc + l)/3i/3fc_i, (52) 

the solution of which is 

/3fc~(fc + l)!/3f (53) 

with /3i = 2. Nicely enough this is the same type of factorial behavior we expect in quantum 
field theory. Moreover, since the coefficient of the factorial is the one-loop /3 function it is 
natural to identify the origin of this behavior with a renormalon. 

In order to clarify the meaning of this renormalon let us consider the perturbative solution 
of equation (HSj) assuming that a <^ 1. Substituting 

/3 = 2«2 (1 + e) , (54) 

in fHS]) we find that e (a) satisfies the equation 

ds s 
2a2— = 3a{l + e) = e-3a + (e^ ea) . (55) 

(JjCx -L ~|~ c 

Because both e and a are much less than unity we can neglect nonlinear terms. Solving the 
resulting linear equation one obtains 

where Pi = 2, Ei (z) is the exponential-integral function and C is the constant of integration. 
Now using the asymptotic expansion of the exponential-integral function at large argument 

-w4(|:^o(^)). (") 

we find that asymptotically the coefficients of the (3 function in the perturbative expansion 
grow as f3k ~ k\(3^. This completely clarifies the origin of the renormalon, which is an 
artifact of the asymptotic expansion of non-analytic function. It also follows from (156|) that 
the accuracy of the Borel resummation does not exceed 

e~W^ ~ — -, (58) 

Re 

which is the expected non-perturbative uncertainty! 
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V. THE NONPERTURBATIVE SOLUTION AND CONFINEMENT 

Finally let us use the exact P function equation to determine what happens in the infrared 
region r > R^, where the perturbation theory is completely out of the control. Obviously 
the nonperturbative effective coupling constant can be used to define a static inter- quark 
potential. Regarding the nature of the classical sources we will leave beyond the scope of 
paper and concentrate mostly on the calculation of the behavior of the running coupling 
constant in the infrared region. As we will see this coupling constant determines-similar to 
QCD- the confining potential between sources. In particular we will compute the energy of 
an isolated source and a dipole, built out of two opposite charges. In the first case we find 
the divergent energy which is an indication of confinement and can be interpreted as the 
absence of isolated sources. The energy of the dipole is finite and positive. Its typical size 
is of order Re and the binding energy ~ R~^ in case of sources with negligible masses. This 
is an encouraging hint towards explaining the colorless hadron states in QCD. 

A. The infrared coupling constant 



Equation (148 p for the exact /3 function can rewritten in the form 



(59) 



da 2a (3 

and can be fully investigated using the phase diagram method. The particular solution we 
need is determined by the perturbative initial condition (!30|) . The resulting nonperturbative 
P function, shown in Fig. [H reaches a maximum value about 0.64 for a ^ 0.69 and after 
that decreases and vanishes at a ~ 0.98. Such behavior of the /3 function in nonperturbative 
regime is quite nontrivial. The most dramatic effect is the absence of any divergence for 
the coupling constant, meaning that in case of exact /3 function, which accounts for all 
resummations, the IR Landau pole is absent. This nonperturbative resummations leads to 
finite running coupling for any finite x interpolating smoothly between the asymptotically 
free UV-regime and the IR- region r ^ Re. 

Beyond the first zero of the /3 function it is more convenient to draw a and /3 separately 
as functions of x = In {r/Rc) because both a and /3 become oscillating functions of the 
scale. In order to find the solution in this region it is more convenient to work directly with 
equation (jH]) instead of ( 159|) . This is the equation for a particle "moving" in a positive 
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FIG. 1: Numerical evaluation of nonperturbative /3 function 



quartic potential in the presence of negative friction. Such particle "oscillates" and if we 
neglect for a moment the friction term, the typical "period of oscillation" can be estimated 
on dimensional grounds as 

Ax - 



where / is the typical amplitude of oscillations. It is clear that this estimate is valid only if 
Ax ^ 1 because otherwise friction term dominates and completely damps the oscillations. 
However, for a (x) ^ 1, when Ax ^ 1, the friction is not so crucial and the system undergoes 
oscillations with the amplitude slowly growing due to this negative friction. In order to find 
how fast this amplitude grows we multiply equation fj4^ by / and rewrite it in the form 

(//')'-/" -^(f)' + «or = 0. (60) 

Averaging this equation over the period of oscillations we find 



Multiplying now f H^ by /' leads to the equation 



(61) 



^r + \aof ] = f\ 



(62) 
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which after averaging and taking into account (16T|) gives us 

d in 4 



dx 



;/*> 



(63) 



Finally solving this equation we obtain 



;r> = cexp(|)=c( 



4/3 



(64) 



leading to the following nonperturbative behavior of the running coupling constant 



a (r) = aof ~ aoA/(/^)cos^ 



0(1) 



2/3 



COS 



2/3 



-Re / V -Re 

for r ^ i?c. In Fig. [2] we summarize the behavior of the running coupling. 

a{r) 



(65) 

As it was 




perturbative 
result 



FIG. 2: Nonperturbative running coupling 

said above although according to perturbation theory the coupling constant should become 
infinite at r = R^ in reality this does not happen. The coupling constant remains finite and 
at scales larger than the confinement scale grows as r^/^. This is the main non-perturbative 
result we can extract from the exact /3 function equation. 
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B. Confinement 

As already discussed we can mimic confinement identifying tlie classical sources as static 
quarks and defining a quench approximation to the static inter- quark potential in terms of 
the effective coupling in the IR region. 

Since in our case the self-interaction contribution to the energy goes like 0^ ~ /^ ~ a{r)'^ 
we can define the static potential as: 

V^(^) ^ ^!M (66) 

r 

Using the non-perturbative value of the running coupling (1^ we get 

V{r)^0{l)R-'(^-^y' (67) 

To check this qualitative result let us compute the energy of the field created by a static 
external charge. Since the field is static and spherically symmetric, the total energy is given 
by the expression 



1 f I /^ i\2 . -^0 ^4 \ „r3^^ r._ f l/c: i\2 . ^0 lA 



E = - I [ (V</))' + ^0^ ]d'x = 2n I ( (a,0)' + -^0^ ) r'dr 



2nQ^J{irdJ-fr-laofy^ 



(68) 



Note that although the contribution of the second term is negative for negative Aq, the total 
energy is positive because the gradient term dominates. The integral above diverges when 
r — 7- 0. This divergence has an ultraviolet origin and it is the same as well known divergence 
of the self-energy of classical point-like electric charge. Therefore it can be removed using 
standard methods. We will focus instead on the IR contribution to the energy. Taking into 
account that at r ^ Re, 

{{rdrff) = {n = ao{n, (69) 

the following expression for the infrared contribution to the total energy is obtained 

^..e^£„„(;.)*.0(l)«^i-(^)"\ po) 

in agreement with the qualitative expectations. The energy of the isolated charge diverges 
as r^/^ and therefore it cannot exist as a free asymptotic state. This can be interpreted as a 
hint of confinement of isolated sources. 
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One can also build ^^ colorless configuration" using two opposite charges Q and —Q sep- 
arated by distance /. At distances r ^ I, the field decreases as r~^ and equation ([3]) 

becomes 

1 d ( ,d(P\ QH^ 

;:i^^-^J=-4-Q^(-) + ^o^' (71) 

Since in this case the anti-screening effect, determined by the last term in this equation, 
is completely irrelevant at large distances we conclude that the total energy of the dipole 
system is infrared convergent. When the distance between charges exceeds the confinement 
scale Re the infrared contribution of the scalar field becomes essential and the total energy 
is 

E-0{l)R-'[^-^ . (72) 

Hence, the interaction potential between two charges grows as distance in power one third. 
This can be interpreted as a confining potential leading to a natural estimate for the mass 
scale of the dipole configuration to be of order m ~ O (1) R~^ . 

VI. DISCUSSION AND SPECULATIONS 

We have shown that certain essential properties of the quantum field theory usually 
considered as having quantum origin can be revealed already at the classical level. In 
particular, the renormalization group structure of the theory including the phenomenon of 
dimensional transmutation is already encoded in the classical equations. 

So far we have considered only the self-interacting scalar field with negative coupling 
constant and external sources. One can naturally ask up to what extent the qualitative 
results obtained in this paper are useful in application to gauge theories, such as QCD. An 
encouraging sign, that classical RG treatment can be generalized for such theories is provided 
by the following simple scaling argument. As we have found, the logarithmic effect of anti- 
screening comes from the term of 0'^ in the equation for the scalar field. In the perturbation 
theory this term represents the density of the charge induced by self-interaction and it drops 
as r~^ as distance r grows. In QCD the gauge field equations for gluons contain two kinds 
of self-interaction terms which drop in a similar way, namely, A^ and ADA. Only ADA 
gives the negative contribution to the /3 function. This term leads to anti-screening effect 
inducing the density of the colored charge decaying as r~^ similar to the case of scalar field. 
Because the structure of the self-interaction terms is different (in one case it is 0^ and in 
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the other A^dA) the interaction potential between two sources in gauge theories can grow 
with the distance not necessarily as r^/^, but as r", where < a < 1. So, the linear growth 
is not excluded. However, the linear growth, although leading to confined charges does not 
necessarily imply the formation of QCD flux tube (see Appendix B). 



Appendix A: On the triviality of Xcp'^ theory with positive A. 

We can use the exact /3 function equation to check the triviality of Ao0^ theory in the 
case of positive Aq in four dimensions (this triviality was rigorously proved in five and higher 
dimensions in J8|). In the case of positive Aq it is convenient to change the signs in the 
definitions of ao and x, so that, 

tto = AoQ^ > 0, X = In (ro/r) . (Al) 



In this case the (3 function defined in (!29|) is related to the standard /3st function used in 
the literature as Pst = /3/2. With these redefinitions the equation for the exact /3 function is 
obtained from (HSj) by substituting /3 — )■ /3 and a — )■ —a : 



2 \da a 
For a ^ 1 the perturbative solution of this equation is 



/3 = 2a2-6a^ + 48a^-570a^ + ... (A3) 



The Gell-Mann-Low equation to one loop, 

da ix) 



dx 
gives us 



2^2 (a;) ^ (A4) 



^' l-2aox l-2Aog2 1n(ro/r)' ^ ' 

According to this result the coupling constant blows up at the Landau pole 

1 



tl = roe 2AoQ^ . (A6) 

The essence of the proof of Xcj)"^ triviality can be reduced to showing that this UV pole will 
survives at nonperturbative level. This is not so obvious because as we have seen the IR 
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one loop Landau pole (for negative Aq) disappears after resummation of the perturbative 
expansion. To find out whether the pole survives or not for positive A, let us solve equation 
(IA2p in strong coupling regime, a ^ 1. We will do it perturbatively in terms of the inverse 
powers of a. Neglecting the linear /3 term in flA2p we have 



^ - — ~ 4a2 (A7) 

da a 

and the corresponding solution of this equation is 

(3 = V2a3/2 (^g^) 

Rewriting ( 1A2I) as 



^ _ ^ = 4«2 _ 2/3, (A9) 

da a 

substituting in the right hand side of this equation the result flASP and solving the obtained 
inhomogeneous linear equation for /3^ we obtain 

(3 = V2a'/' I 1 - ^a-'A = V2a'/' | 1 - ^^ + | ( ^ ) ) ) (AlO) 




This procedure can be repeated recursively giving us higher order power corrections in the 
expansion in Xj \fa ^ 1. Thus we see that for very large a the solution (1A8I) becomes more 
and more accurate and the behavior of /3 functions confirms the expectations in [7|. The 
solution (lASP exactly matches the one loop result in (IA3P at a= Xjl. Using this to fix the 
integration constant in the Gell-Mann-Low equation 

^^ = ^a^l\ (All) 

dx 

we obtain the following non-perturbative result valid for a ^ 1 

« (0 = 8 { ^^ rl ' = 8 { r^^^V^^l ' (A12) 

Thus we see that even after complete resummation of the perturbative expansion the Landau 
pole survives and its non-perturbative location is at (Fig. [3]): 



tl = —e 2Aoe^ . (A13) 

This leads to the triviality of A0^ theory. 
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FIG. 3: Landau pole 
Appendix B: Confinement ^vith and ^vithout strings. 

As we have shown in the paper the potential between the external sources grows un- 
bounded with the separation as a third power of it. In such a picture, charge and anti- 
charge are confined but the confinement in not due to formation of a fiux tube (string-type 
object) but rather due to formation of a finite energy dipole. One may think, that this is a 
peculiarity of not having a linear growth of the potential. In this appendix we will discuss 
this issue and show that even the linear potential does not necessarily imply the existence 
of a string. 

To demonstrate this we begin with an 5*0 (3) sigma model of an isotriplet scalar field 0" 
(a = 1,2,3) with Lagrangian: 



£ = 9^0"9''0" - A' 



,2\2 



(Bl) 



li this case the equation of motion has the static spherically-symmetric solution (see e.g. 

1) 

r = fir)-, (B2) 

r 

where x"' are Cartesian space coordinates and /(r) is the function with the following asymp- 
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totic properties 

/(O) = 0, /(r)|,» (,,)-! ^ V (B3) 

The size of the region where /(r) is different from v is of order re = (Xv)^^. The solution 
above is the 't Hooft - Polyakov monopole in the hmit of zero gauge couphng. Because in 
this hmit the gauge fields become massless and decouple from this solution is often referred 
to as a global magnetic monopole. 

The energy of this monopole can be easily estimated by considering separately the con- 
tribution from core (r < Vc) and from the rest. The core contribution is of order v, and 
can be neglected compared to the energy in the region r > Vc, where /(r) can be set to be 
equal v. Then the contribution of the gradients of the angles (Nambu-Goldstone modes) is 
divergent and we have to cut-off the integral at some R ^ r^. The resulting energy of the 
isolated "charge", 

Er>r. ^ 4:7r drv'^ ~ 47^i^^'^ (B4) 

is linearly divergent. This implies that the potential between two opposite charges is linear. 
In fact, let us consider an anti-monopole placed at distance R from the monopole. The effect 
of anti-monopole is to cut the divergent integral aX r = R and the resulting potential is 

VR^r.^R/rl, (B5) 

thus confining monopole-antimonopole configuration! This picture is very different from the 
QCD flux-tube (string) confinement. To understand this difference let us confront them 
considering heavy quark- ant i- quark pair placed at distance R apart. In the absence of light 
quarks this distance R can be much larger than the QCD scale, tqcd = ^qcd- ■'■^ string 
picture the force between this pair is mediated by a stretched string (electric flux tube) of 
constant tension ~ ^%cdi giving the potential 

yR>rQCD - R/rqcD^ (B6) 

similar to the monopole-anti-monopole potential. However, in the monopole case the flux is 
not confined to a string and for monopole-anti-monopole it has a dipole configuration (see 
Fig. H]). As a result the "monopole color" is not bounded within the string of width Tc and 
can be probed everywhere in the space around the monopole at r > re- On the other hand, 
in the case of string the only possibility to probe the color of the charge is to penetrate 
within the string of width Tqcd- 
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FIG. 4: String vs. dipole 

The difference between these two pictures can be stressed even more if we notice that 
the theory fIBip allows to be "deformed" to the theory in which true strings connecting 
monopoles appear. The appearance of "open-color" monopole-antimonopole configuration 
is due to S{3)/U{1) topology of the vacuum manifold with nontrivial 7r2. To change this 
topology we can further deform the vacuum manifold by spontaneous breaking of the remain- 
ing U{1) symmetry. This can be done by introducing additional scalar field Xa (« = 1, 2) in 
a doublet representation of the SO (3) group. The Lagrangian then becomes 



C = d, 



+ d,x*d^X->^' 



vy - K {X*X - viY + h x*<Px + h'*<P\ , (B7) 



where the contraction of indices is obvious. The parameters {X,Xi,v,Vi,h,h') are chosen 
in such a way that the field x develops an expectation value v' <^ v. In this limit, core 
of the monopole remains nearly unchanged. However at very large distances the field x 
dramatically changes the monopole field. The presence of the second field with nonzero 
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expectation value leads to the following hierarchical symmetry-breaking pattern, 

S0{3)^ f/(l) -^ 1, (B8) 

and the vacuum manifold becomes topologically-trivial. As a result, the static isolated 
monopoles can not exist anymore. However, due the hierarchy of symmetry breaking, v ^ v', 
monopoles do not simply disappear from the spectrum, but rather get connected by the 
strings. Either strings or monopoles are stable in two limits: f = oo and finite v' or v finite 
and v' = 0, respectively. However, when both v and v' are finite, they can only exist as 
hybrid configuration, namely, monopoles connected by strings. This picture is more close 
to the usual string confinement because here the monopole magnetic flux gets confined into 
the string. The thickness of this string is ~ 1/v' ^ Vc, and its tension is v^lnL, where L 
is the string's length. In other words, the tension of the string is logarithmically divergent. 
For example, for the string oriented in z direction the field x near the string but far away 
from the monopoles, is 

Xa ^ Slf (p) e'' , (B9) 

where f{p) vanishes at p = and approaches constant for p > 1/v' in the cylindric coor- 
dinates p, 9. The energy of this configuration is logarithmically divergent with the natural 
cut-off scale of order string size. 

The picture above can be summarized as follows. For R <^ 1/f ', the potential between 
monopoles is linear and field configuration is of dipole type, but for R ^ 1/v' the flux is 
not spread anymore and becomes confined by a string. This leads to the modification of the 
potential f IBSp . which becomes 

V{R)^v'^R\nR. (BIO) 

This is not such a dramatic change in the potential, but more important is the qualitative 
change of the physical picture, because now the monopoles are becoming confined by the 
string. The consideration above illustrates a very important point, that the confinement can 
have very distinct physical origin for the same growing potential between charges. 

Finally, we will consider here one more theory, where confining potential is due to usual 
electric flux. Let us consider U{1) theory with Lagrangian 

L = {F.^F'^T + A,f . (Bll) 
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Unlike (IBip , which describes healthy theory, the legitimacy of this theory as a quantum 
field theory of U (1) gauge field is much less obvious. However, since we are interested only 
in geometric properties of the classical electric fluxes, we will use it for our purposes. In this 
case the equations of motion are 

d>^{2a{FY''F^.) = Ju. (B12) 

For the static charge j^ = 5^6{r)Q^ which produces spherically-symmetric electric fleld, 

F,o ^E,{r) = E{r)^, (B13) 

they become, 

& {2aE{rf''-^Ej{r)) = 6{r)Q. (B14) 

It immediately follows from here that 

2a£;(r)2°-i = ^, (B15) 

and hence 

Q_ 

2ar^ 






Eir) = ( ^ ) . (B16) 



The energy of an isolated charge smeared over a sphere tq diverges as 



R 



Q 



2a 



Echarge= j ^'^^ ( ^ ) '"'^ ~ Q^ ^^ , (B17) 

for either a > 3/2 or a < 1/2 when exponent is positive. Thus, in both these cases 
the energy of an isolated charge diverges. However, in this case the flnite energy of the 
charge-anticharge conflguration is not for granted automatically! The situation is much more 
subtle than in the usual case because of very strong non-linearity, which make superposition 
principle not applicable to Ej. However, the superposition principle in this case is valid for 
E{r)'^°'~''^Ej{r). Therefore, the electric fleld of the dipole of size D is given by. 

The corresponding energy is 

Echarge ~ (Z^Q)^/?^, (B19) 

and it is flnite only for a > 3/2. 
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